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Motivated by the recent re-confirmation by CoGENT of the low-energy excess of events observed 
f^ ' last year, and the recent improved limits from the XENON-100 experiment that are in contention 

r\l , with the CoGENT data, we re-examine the low mass neutralino region of the Minimal Supersym- 

metric Standard Model and of the Next-to-Minimal Supersymmetric Standard Model, both without 
'~{ ' assuming gaugino mass unification. We make several focused scans for each model, determining 

conservative constraints on input parameters. We then determine how these constraints are made 
QQ . increasingly stringent as we re-invoke our experimental constraints involving the dark matter relic 

abundance, collider constraints from LEP and the Tevatron, and then from flavour physics, as a 
series of successive 2a hard cuts. We find that for both models, when all relevant constraints are 
applied in this fashion, we do not generate neutralino LSPs that possess a spin-independent scatter- 
ing cross section in excess of 10~^ pb and a mass 7 GeV < m^ < 9 GeV that is necessary in order 
to explain the CoGENT observations. 
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I. INTRODUCTION 

Recent data from several direct detection (DD) dark matter (DM) experiments have significantly heightened inter- 
ests in the search and identification of the elusive DM, that is predicted to dominate the matter component of the 
Universe, in the form of weakly- interacting massive particles (WIMPs). These include the observations made by the 
CoGENT experiment [1|, ap-type point contact Germanium detector, whose enhanced sensitivity, due to an unprece- 
dented combination of target mass and reduced electronic noise, revealed a possible excess of low energy events [3] 
that has recently been re-afhrmed in June 2011 ^]. If the exponential trend in these low recoil energy events is fully 
attributed to the elastic scattering of WIMP DM, then this signal is consistent with limits on the DM mass in the range 
7 GeV < rriDM S,^ GeV and nucleon scattering cross section in the range 4 x lO""*^ cm~^ < ctjv ^ 1-5 ^ 10~^°cm~^, 
respectively 0, Q ■ Remarkably, the CoGENT results favour a WIMP mass consistent with that indicated by the 
annual modulation count rateQ measured by DAMA/LIBRA [J, Q , assuming elastic scattering and that ctjv is spin- 
independent (SI) fel-fio!]. 

However, the majority of the region in ttidm - cjv parameter space favoured by CoGENT appears to be in cc)ntention 
with recent results from direct detection experiments, in particular, those recently reported by XENON-100 [11], whose 
2011 data consist of non-observations that appear to completely exclude the region favoured by CoGENT [12|. The 
CoGENT and DAMA/LIBRA results are also in contention with the recent observation of two "anomalous" events 
by CDMS-II [i3, ll^) although the discrepancies here are much less extensive than with XENON-100, and such events 
are not statistically significant enough to be able to claim a detection of DM (see, e.g., [15|). Some authors have 
previously concluded that the extent of the discrepancies between the results from CoGENT and DAMA and those 
from the other DD experiments may be partially reconciled by, for example, assuming a different proportion of elastic 
scattering events at DAMA/LIBRA that are channeled, adopting DM halo models (specifically pertaining to the 
escape velocity and velocity distribution of WIMPs) , assuming a lower scintillation efficiency for Xenon, or assuming 
different couplings between WIMPs and either protons or neutrons (see e.g., [ig] and references therein). 

Motivated by past experimental findings, several authors have investigated the parameter space of various particle 
physics models capable of generating light WIMP DM candidates that may potentially explain the CoGENT and 
DAMA/LIBRA observations (see, e.g., |16l - l20| ). However, many such investigations achieve this by conducting focused 
scans with finely-tuned parameters that result in points lying close to the CoGENT/DAMA favoured region, but 
provide a very poor fit to astrophysical and collider constraints. 

In light of the recent release of the 2011 data from CoGENT and XENON-100, we follow a similar line and re- 
examine and directly compare the general behaviour of the neutralino low mass region of two popular particle physics 
models: the Minimal Supersymmetric Standard Model (MSSM) and the Next-to-Minimal Supersymmetric Standard 
Model (NMSSM). We select to investigate the NMSSM in addition to the popular MSSM since in the NMSSM we 
expect that a lower mass lightest supersymmetric particle (LSP) is more easily allowed, relative to that say in the 
MSSM, owing to the weakening of the LEP bounds on the Higgs mass and couplings to Standard Model (SM) particles 
due to its extended Higgs sector. 

To conduct our investigation we use the nested sampling (NS) algorithm, implemented in the Bayesian interface tool 
Multinest [21|, that is incorporated into the SuperBayeS !22| and NMSPEC [23| packages, to perform scans focused 
on the low neutralino mass region when imposing experimental constraints in the likelihood function. However, unlike 
other recent studies, we then also impose the same constraints on the parameter space via a series of {2a) hard cuts 
on the resulting scan data. These hard cuts consist of an increasing number of constraints involving, firstly, those 
on the cosmological DM relic abundance, as inferred from cosmic microwave background (CMB) and galaxy cluster 
observations, secondly, from direct search limits inferred from particle colliders, e.g., LEP and the Tevatron, and, 
thirdly, from indirect constraints from flavour physics and the anomalous magnetic moment of the muon (g — 2)^. 
For each set of hard cuts we highlight the parameter space for each of our models that is consistent with the limits 
on the LSP mass and SI nucleon scattering cross section that are imposed by DD limits. 

The structure of this paper is as follows: In Sec.|ll]we describe the details of the scanning procedure that we utilise 
to investigate the low mass sectors of our two particle physics models. In Sec. lIIII we present and discuss the results of 
our scan of the MSSM obtained when imposing our three sets of constraints from (i) the DM relic density, (ii) collider 
physics (excluding flavour physics) and (iii) indirect limits (flavour physics and {g — 2)^). In Sec. lIVI we present the 
results of our analogous scans on the NMSSM, and compare these results with those corresponding to the MSSM. 
Finally, in Sec.|V]we summarise our conclusions. 



^ In their June 2011 paper, CoGENT also claim to observe evidence of an annual modulation in their data that is consistent with the 
DAMA/LIBRA results ,M- 



II. OUR PROCEDURE 

We are particularly interested in the light neutralino LSP region of a few GeV. Since this window is excluded in 
the usual MSSM by LEP constraints, we relax the assumption of gaugino mass unification. We do the same also in 
our study of the NMSSM, even though in this model a neutralino LSP can potentially evade LEP constraints more 
easily than that in the MSSM, and can therefore be much lighter, because of the possibility of a sizeable singlino 
component. We are particularly interested as to how light the neutralino LSP in the MSSM and the NMSSM can 
be and which experimental constraints dominate in defining the lower limit on its mass. Therefore, in this study we 
explore the low neutralino mass regions of both the MSSM and NMSSM, in each case by performing a single, focused 
NS scan over a wide range of input parameters. 

In TableUwe summarise the sets of experimental and theoretical constraints that we apply to our scans of the MSSM 
and NMSSM. In the top part of TableUwe list the constraints relating to DM relic density and flavour physics that 
we invoke in the respective likelihood functions, using conventional Gaussian probability density functions (pdfs), as 
described in, e.g., [20|. The central part of TableUlists the collider constraints that we include in the MSSM likelihood 
function, again using a Gaussian approximation, associated with our corresponding scan. The exception to this is the 
constraint on the mass of the pseudoscalar Higgs tti^ which takes the form of a non-trivial function firrifn ) (with a 
corresponding inverse function ^(m^)) of the lightest CP-even Higgs mass m^.^ (and tan/?); see, e.g., [2^ for further 
details. We approximate this constraint in our initial scan by invoking the constraint m/^ > 89 GeV in the likelihood 
using a half-Gaussian approximation (we then later invoke the full functional constraints given by /{rrih-^) and g{mA) 
as a 2(7 hard cut; see below). We also invoke all relevant collider constraints in the NMSSM likelihood function 
during our respective scan. These constraints, which are specified in the NMSSMTool code, are applied in our scan 
only as 2(t hard cuts in the relevant likelihood function due to computational challenges. Because of the complicated 
process in which collider constraints are invoked in the NMSSM we omit them in Table|T] and instead refer the reader 
to [2^ for further details. Finally, the bottom part of TableUlists the nuisance parameters, and their corresponding 
ranges of values, over which we perform our scans. Unless otherwise stated, all ranges of values provided in Table|T] 



Category 


Experimental constraint 


Reference 


DM relic density 


n^fe^ = 0.1123 ± 0.0035 (expt) ± 0.1 Q.^K^ (th) 


[24] 


Flavour physics 


BR{B -^ X,7)/10"'' = 3.55 ± 0.26 (expt) ± 0.21 (th) 

BR{Bu -J- TJ/r)/10"* = 1.41 ± 0.43 (expt) ± 0.38(th) 

BR(Bs -> /i+M") < 5.8 X 10"* (95%C.L.) 

5{g - 2)„/10-"' = 29.5 ± 8.8 (expt) ± 1 (th) 


[25] 
[26] 
[26] 

[26] 


Collider 
(MSSM only) 


Finv = 499 ± 1.5 MeV 
m^± > 94 GeV (95%C.L.) 

mjij — g{mA) (2o" hard cut only) 
msj, > 73 GeV (95% C.L.) 
mnj, > 94 GeV (95% C.L.) 
m?, >81.9 GeV (95% C.L.) 

niA = f{mhi) (2o" hard cut only) 


[23 

mi 

[28] 
[27] 
[27] 
[27 
[28 


Collider 
(NMSSM only) 


As implemented in NMSSMTools. 


[23] 


Nuisance 
parameters 


Mt = 171.4 ± 2.1 GeV (expt) 

rribimt)'^'^ = 4.20 ± 0.07 GeV(expt) 

as{Mz)^^'^ = 0.1176 ± 0.002 (expt) 

l/Qem(Mz)*'^ = 127.955 ±0.018 (expt) (MSSM only) 


[27] 
[27] 
[27] 
[27] 



TABLE I. Top: a list of experimental constraints relating to DM relic density and flavour physics that are invoked in the 
respective likelihood functions, using a Gaussian approximation, associated with our scans of the MSSM and NMSSM. Centre: 
a list of collider constraints that are invoked into the likelihood function, using a Gaussian approximation [23] , associated with 
our scan of the MSSM. The exception to this is the constraint on niA and rrih, which we approximate in the likelihood via a 
half-Gaussian corresponding to m^^ > 89 GeV. Then, later, as we do for all other collider constraints, we invoke the constraints 
given by firah-^) and g{mA) as a 2cr hard cut on the points selected by our initial scan. The collider constraints invoked in 
our corresponding scan of the NMSSM are invoked as hard cuts in the respective likelihood function. The procedure in which 
collider constraints are invoked in our scan of the NMSSM via the NMSSMTool code is significantly more complex than that 
which we adopt for the MSSM, and we refer the reader to [23] for further details. Bottom: a list of the nuisance parameters, 
and their corresponding ranges of values, over which our scans are performed. Unless otherwise stated, all ranges of values 
provided are displayed along with their corresponding la experimental (expt) and, where indicated, theoretical (th) errors. 



are displayed along with their corresponding la experimental (expt) and, where indicated, theoretical (th) errors. 
Following our initial NS scans we then successively apply hard cuts to the results using the respective 2a ranges of 
specified combinations of the invoked constraints, where in the case of niA and mhi we invoke the constraints given 
by f{mh-^) and g{mA) as a 2a hard cut. 

An alternative to the above strategy would be to perform three separate scans over both MSSM and NMSSM 
input parameters, omitting subsequent hard cuts, where in each case using a likelihood function invoked with a 
successively increasing number of experimental constraints, i.e., 'likelihood(DM relic density)', 'likelihood(DM relic 
density+coUider)', and 'likelihood(DM relic density+coUider+flavour)'. Despite the efficient nature of such scans, this 
strategy possesses two problems. Firstly, using these three different likelihood functions to investigate the low mass 
neutralino regime can result in posterior pdfs centred upon very different regions of parameter space. For example, 
it is quite possible that the favoured regions resulting from the 'likelihood(DM relic density +collider)' scan will not 
significantly overlap with that from the 'likelihood(DM relic density)' scan, making statistical arguments regarding 
the overall favoured region difficult. Secondly, if we include all experimental constraints into the likelihood function, 
the resulting x^, which, given our assumptions and priors, is given by — 21og[likelihood], can be written as 

xlt^xl + xl + -, (1) 

where xf represents the x^ statistics associated with a particular experimental constraint. In NS scans, Eq. ([1]) allows 
for, and often results in, points being selected which predict values that satisfy some constraints extremely well but 
lie far outside of the usual la range of others, often used in traditional fixed-grid scans, whilst still yielding acceptable 
Xjq^ values. Moreover, such a method lacks transparency as to how the results of its scans correspond with individual 
experimental constraints. 

Since in our scans we apply all experimental constraints as Gaussians in the likelihood functions (with the exception 
being those on ?n^ and m/j; see above), and subsequently make 2a cuts on the corresponding results, we are hence able 
to overcome both of the above issues. However, we emphasise that, because of this we will not provide a statistical 
interpretation of the point distribution of our numerical scans and display all our results as scatter plots, illustrating 
only the density of points, rather than pdfs. We use the NS scanning technique effectively to take advantage of the fact 
that it can perform an efficient scan of a multi-dimensional parameter space in order to obtain a representative sample 
of points, for specified ranges of parameters, allowing us to simply observe the general behaviour of this parameter 
space, rather than to draw statistical conclusions from our scans. 

III. RESULTS IN THE MSSM 

In this section, our goal is to map out the mass range of the light neutralino LSP in the MSSM that is consistent 
with current estimates for the present dark matter relic abundance, as well as constraints, both direct and indirect, 
from collider experiments. It is clear that because existing experimental limits are so stringent, especially those from 
LEP regarding the chargino mass, one needs to consider scenarios that have the potential to evade them in order to 
explore the low mass regime in which we are interested. On the other hand, in order to conduct a numerical scan 
that is manageable, one needs to focus primarily upon the parameters which are relevant to the problem at hand. 

Bearing this in mind, we assume minimal flavour violation and perform a scan, as described in SecHH over the 
following MSSM parameters: 

Ml, M2,//, tan/3, m^,mj-^,TO^^, (2) 

where Mi and M2 are the soft bino and wino masses repectively, fi is the Higgs/higgsino mass parameter, tan/3 is 
the ratio of the vacuum expectation values (VEVs) of the up and down- type neutral Higgs fields, ruA is the mass of 
the pseudoscalar Higgs boson, whilst mj and mj are the left and right hand slepton soft mass parameters. (For a 
more detailed description of the MSSM, we refer the reader to a dedicated review, e.g., [SO]-) We adopt log priors for 
all mass input parameters, and a flat prior for tan/3. The prior ranges of the MSSM parameters over which our scan 
is performed are provided in Tablellll 

We take all left-handed slepton soft mass parameters to be family degenerate (i.e., rriej^ — mji^ — rrir^ — rnj ), and 
likewise for their right-handed partners. On the other hand, gluino and squark masses will not play any significant 
role in determining either the relic density ^^h^ or the spin-independent elastic scattering cross section err , and hence 

we fix them at 1 TeV. Likewise, we fix all the trilinear couplings at the electroweak (EW) scale: At — Aj, — At = A ^ 
—0.5 TeV. In order to focus our scan towards the low mass region of a few GeV, we adopt a log prior distribution of 
the bino mass Mi, which primarily determines the mass, m^, of the lightest neutralino. Also we note that allowing for 
a more generous range of parameters by expanding the lower limits beyond those displayed would result in our scan 



Parameter 


Range 


Parameter 


Range 


bino mass 

wino mass 

/i parameter 

ratio of Higgs doublet VEVs 


0.1 < Ml < 40 

90 < Ma < 150 

90 < /i < 150 

1 < tan /3 < 62 


pseudoscalar mass 

slepton-left mass 

slepton-right mass 


85 < m-A < 600 
70 < mj^ < 3000 
70 < mr^ < 3000 



TABLE II. The prior ranges of input parameters over which we perform our scan of the MSSM. All displayed mass ranges are 
given in GeV. We adopt log priors for all mass input parameters, and a flat prior for tan /3 



finding a disproportionate number of points in the region inconsistent with collider constraints. Note that, because 
we relax the usual assumption of gaugino mass unification (i.e., M3 : M2 : Mi ~ 6 : 2 : 1), and take the prior range 
of Ml well below those of M2 and of n, in order to more efficiently generate light, bino-like neutralino LSP that can 
evade LEP constraints on the invisible Z width and the chargino mass 15|. Another parameter that will play an 
important role in our scan is tan /? since it affects both m^ as well as the fermion-sfermion couplings of the neutralino 
LSP, especially in the low mass regime that we are interested in, and we take a very wide range of its values in our 
scans. 

To start with, in the left panel of Fig. [1] we present the results of our MSSM scan in the m^-J7^/i^ plane. Grey 
points denote all the points resulting from our NS scan, where constraints on fl^h^, collider limits and fiavour physics 
constraints, as specified in TablelH are included in the likelihood. Since we display our results in the m^-n^h^ plane 
we of course omit any hard cuts on the relic abundance, but instead simply illustrate the corresponding 2ct range 
(dashed lines). 
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FIG. 1. Left panel: 



n^h? in the MSSM. Grey points: all the points resulting from our NS scan, where constraints 



on Q,j^h? , collider limits and fiavour physics constraints, as specified in TableU are included in the likelihood. Cyan points: 
points from our scan surviving subsequent 2ct hard cuts using collider bounds. Magenta points: points from our scan surviving 
subsequent 2cr hard cuts using both collider and fiavour physics constraints. We also indicate the 2cr range of relic densities 
within the experimental best-fit value (black dashed curves). Right panel: m^ vs. a^ in the MSSM. Blue points: configurations 
in our initial NS scan yielding relic densities in the range: 0.090 < 0,^^^ < 0.135 (i.e., 2ct range). Green points: points that 
survive a 2a hard cut using both relic density and collider bounds. Red points: MSSM parameter configurations surviving 2a 
hard cuts using all constraints listed in Table|T] We note that the results corresponding to Q-^h'' < Qdm/i^ are very similar 
to those displayed in this panel following our hard cuts, with only a handful of additional points being found in the region: 
m^ > 10 GeV, ap < 2 x 10^^ pb. We also illustrate the regions of parameter space currently favoured by the CoGENT 
(outlined by the magenta dash-dotted curve) and DAMA/LIBRA (with ion channelling, outlined by the dashed grey curve) 
experiments, and display the current limits from CDMS-II (dashed black curve) and XENON-100 (solid black curve). 



Some features are clearly visible. Firstly, we observe a great number of points that predict a relic density that 
far exceeds the corresponding 2a upper limit, particularly those in the low m^ region. However, we also observe 
that there are points in this region that predict a relic density within the 2cr limits indicated by the black dashed 
curves. This is predominantly achieved through efficient i-channel annihilations via light stau exchange, where the 
XT^R coupling, which scales as tan^ /3, and hence, is enhanced at large tan/3. In fact, we find no points with tan/3 < 4 
for m^ < 35 GeV that predict a relic density within the experimental 2cr range. This result is simply a reflection 
of the SUSY neutralino equivalent of the Lee- Weinberg bound on massive neutrinos [3l[. By requiring Vt^K^ to be 
within the 2(7 range, we find ra-^^ > 5 GeV, which is indeed consistent with the results previously obtained in |19| . 
Such low mass points appear to be absent in Fig. 3 of [l^, which finds m-)^ ^ H GeV, although, since cuts were made 
at the 68% C.L., it is possible that points possessing smaller m^ but poorer likelihood were generated there. 

Unfortunately, obtaining points with the desired small values of m^^ of a few GeV comes at the expense of potentially 
grossly violating other constraints, as mentioned in Sec.|TTl In order to see this, we impose a hard cut of 2a on all 
of the collider constraints listed in TablelH The points surviving this cut are displayed in the left panel of Fig.[T] in 
cyan. We can clearly observe that a sizeable number of the points found in our scan (i.e. grey points), especially 
predicting low m^ and/or low ^^h? ^ badly violate collider limits. We also find that, because of the strong correlation 
between the lower limit on m-^ and the minimum slepton mass, invoking the LEP constraint on ?7ie^, m^^ and tti?^ 
are responsible for the observed increase in the lowest value of m^ possessed by a point generated from our scan that 
satisfies the 2a bounds on Vl^h? from roughly 5 GeV (for the grey points in the left panel of Fig.[T]) to approximately 
7 GeV when LEP constraints are invoked (corresponding to the cyan points in the left panel of Fig.[T|). Finally, we 
impose the remaining constraints from flavour physics and the anomalous magnetic moment of the muon, once again 
by making a 2a hard cut on the initial set of points selected in our NS scan. The surviving points are shown in the 
left panel of Fig.[T] in magenta. 

Following the above sets of hard cuts we can clearly observe that whilst the LSP with mass as small as 5 GeV still 
survives, which may lie within the range of WIMP mass favoured by CoGENT and DAMA/LIBRA (with ion chan- 
nelling; see, e.g., f9]), only points with m^ > 13 GeV, which lie outside of the ranges favoured by these experiments, 
predict values of Q.^h? consistent with the corresponding 2a constraint. 

In the right panel of Fig.[T] we demonstrate, this time in the m^ - a^ plane, the effect of applying an increasing 
number of our constraints. The blue points represent the conflgurations in our initial NS scan that yield a value 
of Q,^h'^ within the favoured 2a range 0.090 < il^/i^ < 0.135. The green points illustrate those of the blue points 
that survive a 2a hard cut using collider bounds. Finally, the red points correspond to those MSSM parameter 
conflgurations that survive hard cuts using all of the constraints in TablelH imposed at 2a. 

We note that we also conducted separate scans, identical to those described above except that the neutralino LSP 
generated was allowed to contribute only partially to the currently inferred DM abundance by invoking a 2a upper 
bound of Q.^h? < 0.135 in the likelihood function. The results we found were very similar to those displayed in the right 
panel of Fig.[l] with only a handful of additional points being found in the region: m^ > 10 GeV, a?^ ^ 2 x 10^^ pb. 
The robustness of our results with respect to the constraints on i^-^/i^ is consistent with the findings of Vasquez et 



al. [18|. However we also note that whilst in [18| it was claimed that neutralino mass should obey m^ > 28 GeV in 
order to evade present experimental bounds, determined primarily by the limits from XENON-100 (see Fig. 7 of 18]), 
from Fig.[T] (right panel) we find points with m^ as small as ~ 13 GeV that satisfy both the XENON-100 limits as 
well as our additional 2a hard cuts. 

In the MSSM the SI elastic scattering cross section <7p^ is primarily determined by the i-channel exchange of the 
heavy Higgs scalar H, with the effective xHp coupling being enhanced by large values of tan/3. However, because 
of LEP bounds on the Higgs sector, in the MSSM, one is left with little freedom to boost a?^ towards the regions 
currently favoured by CoGENT and DAMA/LIBRA (with ion channelling), as outlined in Fig. [T] by the magenta 
dash-dotted curve and the dashed grey curve respectively. Moreover, even for the blue points in Fig. [1] (right panel), 
where only a 2a hard cut on fi^/i^ is performed, the Higgs mass is not entirely unrestricted because of the lower 
limits of our selected prior ranges for MSSM input parameters (see Tablelll| . We can also see that imposing collider 
bounds has the substantial effect of cutting off the wedge of points with m^ < 8 GeV possessing the largest SI elastic 
scattering cross sections, a?^ ^ 2 x 10^^ pb, which lie just below the CoGENT favoured region. On the other hand, 
we note that collider constraints have a fairly limited effect at larger values of rriy^ . 

In Fig. [1] (right panel), we observe that when we impose our 2a hard cuts using collider and flavour physics constraints 
a vast number of points in the region: m^ < 13 GeV, (7^^ > 3 x 10~fpb are excluded. Such points are consistent with 
those found in the fore-mentioned investigation of the MSSM in [I^ for 85 GeV < iriA < 100 GeV and tan/3 > 30. 

However, as remarked above, when we invoke flavour physics constraints, particularly those on BR{B — ;■ Xg--^) (and 
also BR{Bs — >■ p'^ii~), see below), as 2a hard cuts, such points are clearly excluded. To confirm this, we note 
that we performed a separate scan of the MSSM, focused on the range of iriA and tan/3 mentioned above and using 
values of the input parameters that range over those utilised in [1^. We found that the resulting points possessed 
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FIG. 2. Left panel: BR{Bs — >■ /i^/i ) vs. a^^ in the MSSM. Grey points: points from our initial NS scan which survive 



subsequent 2a hard cuts using both our Q^^^ ^nd collider constraints. Coloured points: those of the grey points which survive 
2a hard cuts using our flavour physics constraints, excluding that on BR{Bs — > i-i'^jJ.~). The colour scale denotes the value 
of tan/3 associated with each point. We also illustrate the Tevatron limit on BR{Bs — >■ jj.'^fi~) (dashed line). Right panel: 
niA vs. BR{Bs — > I-l'^ fj,~) in the MSSM, using the same colour scale as displayed in the left panel. 



BR{B — i> Xs^) > 4.31 X lO^"* (i.e., 2.3 ct result) and hence, from TableHl fail to survive our 2(t hard cuts, confirming 
the results from our larger scans. 

Upon further examination, we find that it is fairly easy to generate points that satisfy the constraints imposed on 
BR{B — !> Xs7), BR{Bu — >■ tVt), as well as on 5{g — 2)^. In particular, the latter constraint is especially weak since, 
in order to generate a large enough SUSY contribution to the relic density, one requires a fairly light smuon, which, 
being a free parameter, is easy to produce, which in turn can give a substantial contribution to 5{g — 2)^. 

On the other hand, an improved upper bound on BR{Bs ~^ ^^ ^T) from the Tevatron puts a rather strong limit on 
the upper range of a^ values generated in our scan. This can be seen in the left panel of Fig.[21 where we plot the values 
of BR{Bs — !■ fi^iJ,~) against a?^ generated by some of the points in our initial NS scan. Grey points represent those 
scan points which survive subsequent 2a hard cuts using both our DM and collider constraints, while the coloured 
points represent those of the grey points which survive further 2a hard cuts using flavour physics constraints, except for 
that on BR{Bs -^ fj,'^fi~). From Fig.[2]we also observe the well-known correlation between a?^ and BR{Bs -^ iJ,~^fj,~), 
arising from the fact that BR{Bs — > A'^A^") oc tan^ f3/m\ and that uia — mn for ruA ^ 100 GeV [32, |33|- On the 
other hand, at smaller tua (~ 100 GeV) this correlation is relaxed because the approximate mass degeneracy does 
not hold there. 

In Fig. [2] (left panel) we observe two "islands" of points. The one in the bottom left corner is constituted by points 
possessing low to average values of tan/3, which survive the upper limit on BR{Bs — > ii^ ji^) from the Tevatron, 
but produce values of a^^ < 10~^pb, much smaller than those favoured by CoGENT or DAMA/LIBRA (with 
ion channelling). On the other hand, points constituting the upper right-hand island generate values of a^^ up to 
^ 10^^ pb, that include the handful of red points located in the region: 10^^ < a^^ ^ 6 x 10^^ pb, m^ > 13 GeV, 
but are mostly in conflict with the upper limit on BR{Bs -> ^^ jJ.^). Such points possess large values of tan/3 and 
also small values of niA, reaching very close to the corresponding LEP limits for the MSSM, worsening their overall 
likelihood [23, [2^. This can be seen from the right panel of Fig. [51 where we plot the values of tjia against the values 
of BR{Bg -^ ^~^fi~) possessed by the same points as displayed in the left panel. A more thorough treatment of Higgs 
LEP limits by including them in the likelihood could possibly produce more points with a larger a?^ but again they 
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would likely correspond to stretching the LEP limits. 

From Fig. [5] it is clear that in the MSSM, the effect of imposing all the relevant constraints as hard cuts, especially 
that on BR{Bs — ^ 11^1.1^), even at the somewhat relaxed level of 2a, is to exclude any of our points located in either 
of the regions favoured by CoGENT or DAMA/LIBRA (with ion channelling). However, as remarked in Sec. HI we 
note that some authors (see, e.g., [l5| and references therein) have illustrated how these favoured regions may be 
shifted slightly when taking into account certain uncertainties relating to the velocity distribution of Galactic DM or 
the detectors themselves. 

IV. RESULTS IN THE NMSSM 

In this section, we continue our discussion of the low mass neutralino region, now turning to our analysis of the 
NMSSM, starting briefly with a summary of some basic properties of the NMSSM. For a more detailed description of 
the NMSSM we refer the reader to a dedicated review (see, e.g., |34|). 

A. The NMSSM 

The NMSSM is an extension of the MSSM that provides a solution to the so-called /i-problem f35]. Its superpotential 
differs to that of the MSSM in that it contains a new superfield S which is a singlet under the SM gauge group 
SU{3)c X SU{2)l X U(1)y- (We use the same notation for superfields and their respective spin-0 component fields 
for simplicity.) The part of the NMSSM superpotential involving Higgs fields is given by 

W = -e,j\SH:,Hi+^-KS\ (3) 

where Hj ~ {H^,H^), H^ = (7J+,7J°), i,j are SU{2) indices with £12 = —1, while A and k are dimensionless 
couplings in the extended Higgs sector. The superpotential Eq. ([3]) is scale invariant, and the EW scale will only 
appear through the corresponding soft SUSY breaking terms in the soft Lagrangian, £soft • Using the above notational 
conventions, the terms in Csott associated with Eq. ([3|) are given by 

- Csoit = mjf^ H^ Hd + m%^ H^Hu +m^sS*S 
e.jXAxSHlHi + ^kA^ S' + h.c\ , 

(4) 

where rnHj^^Tnu^ and ms are the soft breaking Higgs masses, with A\ and A^ being trilinear soft-breaking terms. 
Consequently, despite that in this model the usual MSSM bilinear /x-term is absent from the superpotential Eq. ([3]) , 
which only possesses dimensionless trilinear couplings, when the scalar component of S acquires a VEV, s — {S), an 
effective interaction ficffH^Hu is generated, with the effective parameter /ieff = As, which is then naturally of the EW 
scale 0. 

In addition to terms from £soft, the tree-level scalar Higgs potential receives the usual D and F term contributions: 

Vf = |Ap (|i?,n5p + \H^m' + \e.,HlH^) + \K\^\St 

-{e,,Xn*WdHiS*^ + h.c.) , (5) 

where gi is the gauge coupling of the electroweak SU(2)i and (72 is the gauge coupling of the U(l)y of the SM. Using 
the minimization conditions derived for the Higgs VEVs, one can re-express the soft breaking Higgs masses in terms 
of A, K, Ax, A^, Vd = (i?°), Vu = {Hi), and s: 



'H. 



A^ (s^ + v^ sin^ /3) - -M| cos 2/3 -I- As tan l3 {ks + Ax) 



^2 _ \2 /„2 , „,2„_2 fl^ , 1 n,r2 



,fj^ - - A^ (s^ + v^ cos^ /3) + -M\ cos 2^ + As cot ^ {ks + Ax) 

2s 



= - A^w^ ~2k^s^ + \Kv'^sm2(3 + ^ sin2^ - kA^s , (6) 



where v"^ = vf + V2 = 2M^/g2, and Mw is the mass of the W-boson. 

Consequently, at the EW scale, the free parameters in the Higgs sector (at tree level) are: A, k, mjj , m?^ , rrig, A\ 
and Ak . Using the fore- mentioned minimization conditions of the Higgs potential, one can eliminate the soft Higgs 
masses in favour of Mz,taTij3 and /i, resulting in A, k, tan/3, /^, Ax, A,^ as a set of independent parameters. In our 
scans we will actually use a somewhat different set of parameters, including the soft scalar masses as well as the soft 
gaugino masses. Mi, M2 and Afs which are all free parameters at the EW scale; these parameters will all be listed in 
Sec. UVBl 

B. Results 

In this subsection we present the results of our scan of the NMSSM which were obtained as described in Sec.lHl 
In Sec. lIIIi we observed that, in the MSSM, even without imposing gaugino unification, the constraints on the Higgs 
sector from LEP and the Tevatron were significant in reducing the number of scan points possessing small values of 
m-^ that also yielded a relic density within 2a of the experimental best-fit value. 

Analogous to our scan over MSSM parameters, in order to evade LEP bounds on the chargino mass for our scan 
of the NMSSM, we again relax the assumption of universality amongst the gaugino masses: Mi, M2 and M3 at 
the unification scale. Consequently, we allow the bino mass to be very small without the implication of yielding an 
unacceptably light chargino. Further, we must also satisfy LEP constraints on the invisible Z width. Because of the 
gaugino non-universality, the majority of the points selected in our MSSM scan were bino-dominated with a gaugino 
fraction gf > 0.6. Similarly, the dominant bino composition of our points selected in our NMSSM scan helps evading 
the constraint on T z, which is proportional to the difference between the squares of the higgsino components of the 
neutralino LSP ^. 

Alternatively, in the NMSSM, one can potentially evade the collider constraints with the neutralino LSP that 
possesses a dominant singlino contribution. Then bino-dominated neutralinos, bino/singlino mixed neutralinos and 
singlino-like neutralinos are possible as a light LSP [Sg. In these three scenarios, for neutralino masses below ~ 
15 GeV, efficient annihilation is achieved mainly via a light Higgs resonance exchange [37|. Therefore, in order to 
obtain an acceptable range of the relic density for a very light neutralino, one needs in general very light Higgses 
with the mass of the lightest Higgs scalar ruh^ tuned to be close to 2m^. This implies that, in order to evade LEP 
bounds, the lightest Higgs must be singlet-like. In order to generate this kind of Higgs boson without driving down 
the mass of the SM Higgs below the LEP bound of ~ 114.4 GeV [2^, one needs to decouple the singlet and doublet 
Higgs components (i.e., obtain the limit -M5i3,A^|23 ~^ ^ for the mass matrix Mg of CP-even Higgs states (see, 
e.g., Eq. (2.8) of [Ml). 

As discussed in [3^, there are two possible ways to achieve this. One is to adopt the limit A — > 0. Alternatively, 
one can enforce the relation: 

2 

^A = . ^o A^cff - 2ns. (7) 

sm 2p 

Adopting Eq. ([7]) exactly yields a mass for the lightest CP-even singlet Higgs scalar given by 

m? =A2^a^ + ks(A, + 4ks), (8) 

^ As 

which is equal to M'g 33 (once again, see, e.g., Eq. (2.8) of [3^). Such solutions were found in the study conducted by 
Belanger et al. [331 ■ From Eq. ([8]) , we can see that there are new types of solutions present in the NMSSM that are 
not obtainable in the MSSM. In the limit A — )► 0, the mass of the CP-odd neutral scalar Higgs (i.e., the pseudoscalar 
Higgs) is given by 



m 



2 =\(2k + —''-^ 

pseudoscalar y\y^rv -r 



= A(2k + — )-u^ sin 2^ - 3kA^s -^ -3kA«s, (9) 



which is equal to 7Wp22' where A4p is the mass matrix of CP-odd Higgs states, and where, in the limit A — > 0, 
M. PA-2 — ?> and hence, due to the zero mass terms of the Goldstone boson, M.-p is diagonal (see, e.g., Eq. (2.10) 
of (S^). From Eq. © we can see that a negative value of kA^s is necessary in order to obtain "T-pgeudoscaiar > ^^'i' 
consequently, avoid tachyonsQ For small A it will be difficult to achieve a resonance annihilation via a singlet Higgs 



^ Notice that without lost of generality A can be chosen to be positive [3S|| and that with this convention the sign of s and the sign of /^cft 
are the same. 
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Parameter 


Range 


Parameter 


Range 


bino mass 


0.1 < Ml < 30 


CP-odd neutral scalar Higgs mass 


85 < mA < 600 


wino mass 


90 < Ma < 500 


slepton-left mass 


70 < mj^ < 3000 


/^cff parameter 


90 < /^cff < 500 


slepton-right mass 


70 < mr < 3000 


ratio of Higgs doublet s 


2 < tan/3 < 65 


trilinear terms 


\A^\ < 100, \A\ < 4000 


Higgs sector coupling 


10"'' < A < 0.5 


Higgs sector coupling 


10"'' < K < 0.5 



TABLE HI. The prior ranges of input parameters over which we perform our scan of the NMSSM. All displayed mass ranges 
are given in GeV. We adopt log priors for all input parameters except tan/3, for which we use a fiat prior. 

when the neutralino is singlino-like due to tendency for the masses of the neutralino LSP and the Higgs to be similar. 
Hence, in this case, a bino-dominated neutralino LSP will be favoured. 

Conversely, for larger A, it is possible to generate nearly pure singlet-like Higgs states by fine-tuning A\ according 
to Eq. ([7]) . These solutions allow for singlino-like neutralinos that annihilate via a singlet-like Higgs that are much 
more efficient than those corresponding annihilations involving a bino-like LSP. Since the constraints from LEP on 
the chargino mass ultimately yields a lower bound on /ioff — As, and from Eq. ([7]), a hierarchy is established between A 
and K, where k <^ X for A 3> 0, in such a way as to result in a small value of m^ '^ 2ks. Alternatively, for large A, one 
may avoid this hierarchy as well as generate a light, pure, singlet-likc CP-even or CP-odd Higgs by finely-tuning A,^ to 
drive down the mass given by Eq. ([5]) or Eq. ^ respectively. Such fine-tuning allows for LSPs similar to those found 
in [20|. As a result of these two fine-tuning scenarios, the solutions in the regime A > 0.1 are generally finely-tuned, 
resulting in mixed bino-singlino LSPs. 

Taking into account all our constraints (i.e., DM relic density, collider and fiavour physics), one can determine that 
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FIG. 3. Left panel: m^ 



in the NMSSM. Cyan points: all the points resulting from our NS scan, where constraints 



on Q^h^ and fiavour physics bounds are included in the likelihood as Gaussians, and collider constraints are included in 
the likelihood as hard cuts. All relevant constraints are specified in TablelH Magenta points: points from our initial scan 
surviving subsequent 2a hard cuts using flavour physics constraints. We also indicate the 2a range of relic densities within 
the experimental best-fit value (black dashed curves). Right panel: rux vs. a^^ in the NMSSM. Green points: those of the 
cyan points in the left panel which survive a subsequent 2a hard cut on Q-^h . Red points: those of the magenta points in 
the left panel that survive a subsequent 2a hard cut on Q^h^. We also illustrate the regions of parameter space currently 
favoured by CoGENT (outlined by the magenta dash-dotted curve), DAMA/LIBRA (with ion channelling, outlined by the 
dashed grey curve), and illustrate the current limits from CDMS-II (dashed black curve) and XENON- 100 (solid black curve). 
We also display the estimated upper bound on ap in the NMSSM according to Das & EUwanger 'iTj] when using default values 
of strange quark content in nucleons (orange curve). 



11 




NMSSM, u „>0 

' ett 

0.090 < a h^< 0.135 

X 



20' 



>^^- 



20 40 60 80 100 120 140 



m^ (GeV) 







* . 



NMSSM, u „>0 
ett 

0.090 <0 h^< 0.135 
X 



20 40 60 80 100 120 140 

m, (GeV) 



h 



N^ 



0.2 



0.4 



0.6 



0.8 



0.2 



0.4 



0.6 



0.8 



FIG. 4. Left panel: mui vs. m^ in the NMSSM, for the points from our initial NS scan that survive subsequent 2a hard cuts 
on Q-^h^ and when using constraints from colhders (including flavour physics). The colour scale illustrates the singlet fraction 
Sg of the lightest Higgs. Right panel: Identical to the left panel except now the colour scale indicates the singlino fraction of 
each respective neutralino LSP. 



a value of A ^ 10~^ favours a bino-doniinated neutralino LSP. For larger values of A ~ 10~^ neutralinos are a mixture 
of both bino and singlino. For this reason, the selection of either a log or flat prior for A can make a huge difference 
to numerical results. We highlight the caveat that all such solutions are fine-tuned since 2m-)^ « to/u or 2m^ « rua^ , 
where ruai i^ the mass of the lightest doublet-like pseudoscalar Higgs (whereas the heavier doublet-like pseudoscalar 
mass is denoted by Waj), however for A > 10~^ we find that a further numerical fine-tuning, this time involving A\, 
is necessary in order to generate a singlet Higgs pure enough to escape collider constraints. One should also mention 
that in this scenario collider constraints on the heavier, MSSM-like, CP-even Higgs /i2 become extremely important, 
favouring small values for tan /? ~ 3 |37i] . 

Bearing in mind the above, we performed our scan of the NMSSM, using the NMSPEC code [23| and imposing all 
relevant constraints in the likelihood as described in SecHH taken over the following input parameters: 



Ml , M2, /icff, TOA, mj^ , mj^ , ^K, ^, tan /3, A, k. 



(10) 



where ttia denotes the running mass of the lighter CP-odd Higgs-doublet, as defined in the NMSPEC code [23| (see 
Table|TTT|, and the other parameters in Eq. ([T0|) have already been defined abovqfl- where niA is an approximative value 
for TTiai- Once again, we adopt log priors for all input parameters except tan/3, for which we use a flat prior, when 
conducting our scans. In particular, we utilise a log prior for A and k so that we may explore the general behaviour of 
low mass LSPs in the A-k plane, rather than adopt flat priors and focus on flnely-tuned solutions possessing A ^ 0, 
as in, e.g., [37]. Also, for the same reasons as, and in order to be consistent, with our scan of the MSSM, we assume all 
slepton soft mass parameters as well as their right-handed partners to be degenerate and flx the gluino and all squark 
masses at 1 TeV. The prior ranges of the NMSSM parameters Eq. pUj) over which our scan is performed are provided 
in Table lHll In analogy to our MSSM scan, we restrict the bino mass to the range 0.1 GeV < A/i < 30 GeV and the 



Note that in the NMSPEC code, when m^ is chosen as input. Ax is treated as output. 
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FIG. 5. Left -panel: m^ vs. a^^ in the NMSSM, for the points from our initial NS scan that survive subsequent 2a hard cuts 
on Q^h^ and when using constraints from coUiders (including flavour physics). The colour scale indicates the singlino fraction 
of each respective neutralino LSP. The experimental bounds and theoretical limits are identical to those displayed in the right 
panel of Fig.O Right panel: niai vs. m^ in the NMSSM, using the same points as those displayed in the left panel. 



higgsino mass parameter to the range 90 GeV < /icff < 150 GeV, in order to focus on light, bino-like neutralino LSP 
and evade collider constraints on the Higgs sector, as discussed above. 

We present the results of our scan in the left panel of Fig. [31 displayed in the m^-Q,^h? plane. Cyan points 
correspond to those resulting from our initial scan, where collider constraints have been invoked in the relevant 
likelihood as hard cuts, differing from our MSSM scan, as described in Sec.|TTl but with the constraints on Vt^h? and 
those from flavour physics being invoked using a Gaussian assumption, as with our MSSM scan. The magenta points 
in the left-panel of Fig. [3] then correspond to those of the cyan points that survive a subsequent 2cr hard cut utilising 
flavour physics constraints. For both data sets, we observe that it is possible to generate points in the NMSSM 
that yield a relic density within 2a of the experimental central value from WMAP, which are indicated by horizontal 
black dashed lines. However, as in the MSSM case, we can see many points where their associated relic density is 
many orders of magnitude larger than 0.1, due to premature thermal freeze-out brought about by the insufficient 
annihilation rates of the respective LSP cases. 

In the right panel of Fig.[3l we demonstrate the effect of applying our different hard cuts on the m^ - a?^ plane. 
The light green and red points correspond to those points that survive respective 2cr hard cuts on Q,^h? on the cyan 
and magenta points displayed in the left panel of Fig. [3] For comparison, we also superimpose current limits from 
CoGENT, DAMA/LIBRA (with ion channelling), CDMS-II and XENON-100. We observe that the majority of our 
selected points typically generate values of the SI elastic scattering cross section in the range a^ < 10^^ pb. We also 
note that many of the green points generate much larger values of Up^ up to 2 x 10"^ pb, however, the bino fraction 
of these neutralino LSPs are too small to survive our hard cuts involving flavour physics constraints. 

From Fig. [3] (right panel) we also observe that our results, following our 2(T hard cuts, coincide well with the 
corresponding upper limits on a^ estimated by Das & Ellwanger [17] (orange curve) using their default model regarding 
the strange quark content of nucleons. Despite this agreement, we also note that the studies conducted by Vasquez 
et al. 18] and Cao et al. [20] generate results that include a small population of neutralino LSPs with a^^ as large 
as 10~^ pb with m^ ^ 10 GeV. We note that the a?,^ values of these points grossly violate the DD limits from 
XENON-100 by up to four orders of magnitude, as do the majority of the points we obtain possessing a?^ > 10^^ pb. 
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FIG. 6. Lejt panel: BR(Bs — >■ fi'^/i') vs. a^ in the NMSSM. Grey points: points resulting from our initial NS scan that 
survive a subsequent 2a hard cut on Q,xh . Coloured points: those of the grey points displayed which survive further 2a hard 
cuts from colliders (including flavour physics constraints, except that relating to BR{Bs — > fj,'^ jj.~)). The colour scale displayed 
indicates the value of tan/3 associated with each neutralino LSP. Right panel: nia^ vs. BR{Bs — J- fi'^ fj,~) in the NMSSM, using 
the same points as displayed in the left panel. 



and hence their omission do not effect our conclusions regarding the viable regions of NMSSM parameter space that 
are consistent with such current experimental limits. Moreover, we also remind the reader that, as described in SecHIl 
we implicitely invoke our collider constraints as 2a hard cuts during our scan of the NMSSM, and hence omit points 
that grossly violate such constraints. 

In Fig.m we plot the respective values of mh^ against m^, associated with the points from our initial scan surviving 
subsequent 2a hard cuts using flavour physics constraints. We observe that, for m^ < 15 GeV, the reason why these 
points survive our hard cuts on fl^h^ is predominantly because of the existence of an s-channel resonance annihilation 
of the neutralino LSP via the lightest singlet-like Higgs, i.e., where 2m-^ ~ ruhi- This can be seen in the left panel 
of Fig.m where we plot the singlet fraction, 5*1, of the lightest Higgs as a third axis. The points corresponding to 
the LSP annihilating through the Higgs resonance also possess the lightest Higgs that is mostly singlet, which is 
essential in order to evade collider constraints on mh^ . We also observe points with m^ > 15 GeV that also survive 
our hard cuts on ft^h'^. These points correspond to those configurations where one has very light slepton masses for 
all three generations, by virtue of our relaxation of the universality of slepton masses, allowing one to obtain a small 
enough relic density through additional t-channel annihilations, as is the case in the MSSM. For masses of the lightest 
Higgs in the range rrihi •^ 114.4 GeV we also observe the onset of the MSSM-like scenario, where the lightest Higgs 
is completely doublet-like, since we are now in the LEP-allowed region. 

In Fig.[S]we display our results in the m^ - Cp plane (left panel) and in the rua^ - m^ plane (right panel). In the left 
panel of Fig.[51 we can clearly see that there are only a minority of our selected points that possess a significant singlino 
component that survive our cuts and generate a value of cr^ > 10^^ pb. In the right panel of Fig.[5]we observe that for 
some surviving configurations the lightest pseudoscalar is lighter than the lightest neutralino. For these configurations, 
if the neutralino LSP also possesses a significant singlino component, annihilations where pseudoscalars are in the final 
state become significant. This explains the presence of those few singlino-dominated points in the right panel of Fig. 21 
where we present the singlino fraction, N"^, of the neutralino LSP of the points displayed in the left panel of Fig.|4] 
(via a colour scale), residing below those undergoing the resonance annihilations via the lightest scalar Higgs (i.e., 
2m^ ~ nihi). Alternatively, there is the possibility of configurations possessing a neutralino LSP that undergoes a 
resonance annihilation via the pseudoscalar Higgs (i.e., 2m^ ~ mai), which is also evident in Fig. [5] Such annihilation 
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FIG. 7. Left panel: logio A vs. logio k in the NMSSM, for the points from our initial NS scan that survive subsequent 2cr hard 
cuts on flxh^ and when using constraints from colhders (including flavour physics). The colour scale illustrates the singlet 
fraction S| of the lightest Higgs. Right panel: Identical to the left panel except now the colour scale indicates the singlino 
fraction of each respective neutralino LSP. 



channels are evidently crucial in yielding an overall annihilation rate sufficiently large enough to generate the correct 
relic density (within the 2cr range) necessary to survive the invoked hard cuts. 

The consequences of imposing the constraints from flavour physics are much more severe. This can be seen in 
Fig.ini where the grey points represent those selected by our initial NS scan that survive a subsequent 2a hard cut 
on fi^ft,^, and the coloured points represent those of the grey points that survive further 2a hard cuts using flavour 
physics constraints except for that relating to BR{Bs -^ ^'^fj,~). In the left panel we display these points in the 
BR{Bs -^ ^~^ fi~) - a?^ plane, and, in the right panel, the Waa -BR[Bs — >■ jX^fjC) plane. In both plots we display the 
value of tan/3 associated with each configuration via a colour scale. 

We observe that, like in the MSSM (see left panel of Fig. [2]), configurations with large tan/3 typically yield a large 
SI elastic scattering cross section within the range 10~^pb < a?^ < 10~^pb, but are excluded by the constraint on 
BR{Bs -> n'^fi^) from the Tevatron. In the right panel, we also observe a correlation between those points possessing 
large values of tan/3 and small values of ttiq^ < 300 GeV and those possessing large values of BR{Bs -^ fi~^iJ,~), 
analogous to the correlation between BR{Bs -^ /i+/x^), niA and tan/3 in the MSSM displayed in Fig. [21 However, in 
the NMSSM, this correlation is evidently much weaker than in the MSSM due to the additional degrees of freedom 
relating to the two pseudoscalar Higgs masses. 

Finally, in Fig. [7] we present our results in the X- k plane with third axes displaying the singlino composition of the 
lightest Higgs scalar (left panel) and singlet composition of the lightest Higgs (right panel) . We clearly observe that 
our scan points are segregated into several different regions in the X- k plane, each with a distinctive minimum value 
of logio -^j containing points that we can ascribe to one of approximately four different categories that we will now 
discuss. 

Starting at small values of A we observe the following categories of points: (i) At A ~ 10^** we can see the onset of 
the MSSM-like region, consisting of bino-dominated LSPs with m^ > 7 GeV (see left panel of Fig. S]), with the lightest 
Higgs scalar being doublet-like with mhi ^ 114.4 GeV in order to evade LEP bounds. In this case the neutralino 
annihilates mainly through a i-channel slepton exchange as described earlier in this section, (ii) For A ~ 10~^ and 
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K ~ 10"'' we observe the onset of a region where m,^ > 15 GeV and, as with the MSSM-hke region, the LSP is 
bino-dominated and stiU annihilates mainly via i-channel slepton exchange, but in this case the lightest Higgs is 
singlet-like, (iii) For A ^ 10"^ we find the onset of those points involving bino-dominated neutralino LSPs that 
annihilate through the fore-mentioned s-channel resonance involving a singlet-like Higgs, possessing a mass following 
the relation: m^ « 2771^^. iv) Lastly, for A > 10~^ we find neutralinos with m^ < 15 GeV which are now a mixture 
of both singlino and bino and, as with (iii), annihilate via a resonance involving a singlet-like Higgs. As explained in 
Sec. lIVBl a degree of fine-tuning in A\, via Eq. (O, or in ^„, via Eq.®, is necessary in order to generate very light 
singlet-like Higgses for large values of A ~ 1. Such solutions can potentially generate large values of the SI elastic 
scattering cross section that are several orders of magnitude larger than that associated with the bino-dominated case, 
for A ~ 10~^, depending on the degree of fine-tuning involved. However, as remarked earlier in this section, such 
points are clearly excluded by the DD limits, particularly those from XENON-100. 



V. SUMMARY 

In this paper we have re-examined the low mass neutralino region of both MSSM and NMSSM parameter space 
in light of the recent re-confirmation of a low energy excess of events observed by CoGENT, and the recent limits 
arising from the non-observations by XENON-100 that are in contention with the CoGENT data. To conduct our 
investigation we utilised focused NS scans when imposing a variety of experimental constraints, associated with (i) the 
cosmological DM relic abundance, (ii) direct collider search limits and (ii) collider physics relating to flavour physics 
and S{g — 2)^, in the likelihood. Unlike the majority of previous studies we then re-invoked each of these constraints 
on the resulting scan data, this time, via a successive series of 2a hard cuts. 

Firstly, regarding the MSSM, we relaxed the unification of gaugino masses in order to efficiently generate light, 
bino-like neutralino LSPs without violating collider constraints involving the invisible Z width and the chargino mass. 
We observed that for many configurations from our initial scan producing the neutralino LSP mass m^ > 5 GeV 
the corresponding values of il^h^ were small enough to be within the WMAP 2a range. This could be achieved 
predominantly through efficient t-channel annihilations via stau exchange, which are enhanced for points possessing 
large tan/3. However, owing to LEP constraints on the slepton masses, when we invoked collider constraints as 2a 
hard cuts, only bino-like LSPs with tti^^ ^ H GeV survived. Moreover, when invoking the constraints from flavour 
physics, the minimum LSP mass was further increased to approximately m^ > 13 GeV, which is slightly heavier 
than the values within the favoured CoGENT region, and much lighter than the 28 GeV limit claimed by Vasquez 
et al. [l^ . We find that for these points it is extremely difficult to generate values of the SI elastic scattering cross 
section much larger than 10~^ pb, owing primarily to the constraints on the x^p coupling, which is enhanced for 
large tan/3, coming from flavour physics, particularly the limit on BR{Bs -^ ^^ij.~). We emphasise that this limit 
is particularly stringent within the framework of minimal flavour violation which we have adopted but can be fairly 
easily relaxed when one abandons it; see [40, |4l|. However, we point out that solutions possessing such large values of 
(Tp\ within our permitted mass range of m^ > 13 GeV, are also clearly excluded by the DD limits from XENON-100. 

Regarding the NMSSM, once again, we conducted a focused NS scan in the low LSP mass region, where (non- 
flavour) collider constraints were invoked in the likelihood this time as 2a hard cuts. Once again, in order to efficiently 
generate light neutralino LSPs without violating LEP bounds on the invisible Z width and the chargino mass, gaugino 
unification was relaxed. Whilst the majority of points surviving collider constraints were found to be bino-like, we 
found that a significant proportion of these points also evaded collider constraints by adopting a partial singlino 
component, made possible through the extended Higgs sector of the NMSSM. In such cases, we found that for light 
LSPs possessing m^ < 15 GeV we were able to generate values of 17^ /i^ small enough to lie within the WMAP 2a 
range through s-channel resonance annihilations via the lightest CP-even Higgs. Consequently, in order to evade LEP 
bounds, in such cases the lightest Higgs had to be singlet-like, very much decoupled from the doublet Higgs. This 
could be achieved by adopting the limit A — >■ 0, in which case bino-dominated LSPs are favoured, or by finely-tuning 
A\ , generally resulting in a combination of bino-singlino LSP configurations and a hierarchy between A and k where 
K <C A for A ::^ that we explicitly described, or by finely-tuning A^^. However, for m^ > 15 GeV, we found that 
points can generate il^h^ small enough to survive our cuts by virtue of additional i-channel annihilations via very 
light sleptons, possible through our relaxation of the universality of slepton masses, or in some cases through s-channel 
annihilations via the pseudoscalar Higgs. We found that the points from our initial scan generated values of ap^ up to 
approximately 2 x 10~^ pb. However, we found that many of these points possessed too large a bino fraction to survive 
the fiavour physics constraints, particularly that coming from BR{Bs — > fi'^ii~), that we later re-invoked as 2a hard 
cuts, finding results very similar to those obtained in [I4I where Cp' < 10^^ pb. We also highlight that many of these 
points with m-^ > 10 GeV are excluded by DD limits from XENON-100, like many of the points discussed by other 
authors, including [l^, l20|, that possess 17^' as large as 10~^ pb, but are generated without invoking experimental 
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constraints as 2cr hard cuts, as we do, in order to omit points that grossly violate individual constraints. 
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